Constraining viewing geometries of pulsars with 
single-peaked 7-ray profiles using a 
mult i wavelength approach 

A S Seyffert 1 , C Venter 1 , T J Johnson 2,3 and A K Harding 2 

1 Centre for Space Research, North- West University, Potchefstroom Campus, Private Bag 
X6001, Potchefstroom 2520, South Africa 

2 Astrophysics Science Division, NASA Goddard Space Flight Center, Greenbelt, MD 20771, 
USA 

3 Department of Physics, University of Maryland, College Park, MD 20742, USA 
E-mail: 20126999@nwu.ac.za 

Abstract. Since the launch of the Large Area Telescope (LAT) on board the Fermi spacecraft 
in June 2008, the number of observed 7-ray pulsars has increased dramatically. A large number 
of these are also observed at radio frequencies. Constraints on the viewing geometries of 5 
of 6 7-ray pulsars exhibiting single-peaked 7-ray profiles were derived using high-quality radio 
polarization data Q]. We obtain independent constraints on the viewing geometries of 6 by 
using a geometric emission code to model the Fermi LAT and radio light curves (LCs) . We find 
fits for the magnetic inclination and observer angles by searching the solution space by eye. Our 
results are generally consistent with those previously obtained [I], although we do find small 
differences in some cases. We will indicate how the 7-ray and radio pulse shapes as well as their 
relative phase lags lead to constraints in the solution space. Values for the flux correction factor 
(fn) corresponding to the fits are also derived (with errors). 



1. Introduction 

We study 6 pulsars (J0631+1036, J0659+1414, J0742-2822, J1420-6048, J1509-5850, and 
J1718— 3825) detected by Fermi LAT, all exhibiting single-peaked pulse profiles within current 
statistics at energies > 0.1 GeV. These pulsars are also detected in the radio band. Both of these 
properties aid in constraining the possible solution space for the respective pulsar geometries, 
as derived from predictions of their light curves (LCs). 

A study of the pulsars' geometric parameters, the inclination and observer angles a and £, 
has been performed [I] using the radio polarization and LC data. They constrained the solution 
space for a and j3 of 5 of the individual pulsars, with the impact angle (3 = ( — a, using fits to 
these radio data as well as predictions for the value of the half opening angle, p, of the radio 
beam derived from the radio pulse width (e.g., [2]). 

The aim of this study is to obtain similar constraints on a and £ of the same pulsars using 
an independent, multiwavelength approach. We use a geometric pulsar emission code (e.g., [3]) 
to model the Fermi LAT and radio LCs, and fit the predicted radio and 7-ray profiles to the 
data concurrently, thereby significantly constraining a and This also allows us to test the 
consistency of the various approaches used to infer the pulsar geometry. 



2. Model 

We use an idealized picture of the pulsar system, wherein the magnetic field has a retarded 
dipole structure [1] and the 7-ray emission originates in regions of the magnetosphere (referred 
to as 'gaps') where the local charge density is sufficiently lower than the Goldreich-Julian charge 
density [5]. These gaps facilitate particle acceleration and radiation. We assume that there are 
constant-emissivity emission layers embedded within the gaps in the pulsar's corotating frame. 
The location and geometry of these emission layers determine the shape of the 7-ray LCs, and 
there exist multiple models for the geometry of the magnetosphere describing different possible 
gap configurations. 

We included two models for the 7-ray emission regions in this study, namely the Outer 
Gap (OG) and Two-Pole Caustic (TPC) models. In both the OG and TPC models (the Slot 
Gap model [6] may serve as physical basis for the latter), emission is produced by accelerated 
charged particles moving within narrow gaps along the last open magnetic field lines. In the 
OG model [7J, radiation originates above the null charge surface (where the Goldreich-Julian 
charge density changes sign) and interior to the last open magnetic field lines. The TPC gap 
starts at the stellar surface and extends along the last open field lines up to near the light 
cylindeiQ where the corotation speed approaches the speed of light [8] . The special relativistic 
effects of aberration and time-of-flight delays, which become important in regions far from the 
stellar surface (especially near the light cylinder), together with the curvature of the magnetic 
field lines, cause the radiation to form caustics (accumulated emission in narrow phase bands). 
These caustics are detected as peaks in the observed 7-ray LC [8j [9] . 

We used an empirical radio cone model pl)|, where the beam diameter, width, and altitude 
are functions of the pulsar period P, its time derivative P, and the frequency of observation v, 
in order to obtain predictions for the radio LCs. This is different from the approach taken by 
PQ, as we have a different prescription for the radio emission altitude and p, and we do not use 
the polarization data. 



3. Method 

Using the above geometric models, we generated LCs as a function of a and (, keeping the gap 
widths and extents constant. Due to the large size of the (a,£)-space at a 1° resolution, it is 
impractical to search blindly for the optimum LC by eye. We therefore started by generating a 
so-called atlas of LCs with a 10° resolution for each pulsar and first identified candidate LC fits. 
We then performed refined searches in the regions of (a,£)-space around those candidate fits at 
a 5° and later 1° resolution. This approach of generating all the LCs in a candidate region and 
comparing them directly, allowed us to obtain constraints on a and We then inferred contours 
for a and C for each pulsar and for each geometric model. We lastly used these contours to find 
the corresponding values of fa (with errors; see Section 3.2) predicted by the OG and TPC 
models. 



3.1. Obtaining the contours 

To obtain an (a,£)-contour delineating the region of plausible LC fits, we start at the candidate 
solution (5° resolution) and first move away from it in steps of 1° in a. At each step we then 
step through C, until the upper and lower bounds of the contour are located for the fixed a. This 
procedure is repeated until no further acceptable LCs are found. As an example, the red regions 
A and B in the left panel of Figure [T] indicate the resulting contours obtained in this manner 
for PSR J1509-5850 for the TPC case. 

1 The slight difference in transverse polar position of the gaps with respect to the magnetic axis (as motivated by 
physical models) will result in a systematic shift in best-fit (a,£)-contours between the two models. See SectionBl 
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Figure 1. Left: The (a,£)-contours (regions A and B) for radio and 7-ray LC fits of PSR 
J1509— 5850 for the TPC model (red). The dark grey and light grey regions correspond to 
double-peaked and single-peaked radio LCs respectively. The solid lines indicate approximately 
the locations of the boundaries to the contours due to: (a) the relative phase lag between the 
radio and 7-ray peaks being too large, (b) the bridge emission's relative intensity being too 
low, and (c,d) the 7-ray peak being too narrow. The dashed line (bounded area) indicates the 
resulting errors on a and £, with the green block indicating our best fit. Right: Representative 
LC fits corresponding to the contours A and B. 



Oftentimes we find two disjoint solution contours, with one of the two sometimes producing 
a better LC fit. The example shown here is one of these cases, as can be seen from the 
representative LCs shown in Figure[TJ Contour B yields a considerably better fit than contour A, 
and thus we can ignore contour A in favour of contour B as indicated by the dotted line in 
Figure [T] The values of a and £ reported in Table [T] are those at the centre of this dotted box, 
while the errors on these values are conservatively chosen to encompass the dimensions of the 
box. In this case the resulting values are a = (61° ± 5°) and £ = (44° ± 7°). 

3.2. Finding /n 

It is important to be able to convert the observed energy flux of a pulsar to its all-sky luminosity. 
The flux correction factor is used for this purpose. It is a highly model-dependent parameter, 
and allows us to determine what fraction of the all-sky luminosity we are observing from a pulsar 
if the geometry is known (a and £), for a given model. It is also crucial in deriving the 7-ray 
efficiency of a pulsar, and is defined as [11] 



fjF 7 (a,C,0)sinCdC# n s 
2 J F 1 {a,CE,(f>)d(f> 

with F-y being the photon flux per solid angle ('intensity'), and Ce the Earth line-of-sight. The 
value of fa is typically taken to be 1, meaning that the observed energy flux is assumed to be 
equal to the average energy flux over the entire sky. Due to the OG model predicting emission 
over a relatively smaller region of phase space, its associated fn values are typically smaller 
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Figure 2. Contour plot of /q as function of a and £ implied by the TPC model for PSR 
J1509-5850. The red area indicates the best-fit (a, C)-contour for the radio and 7-ray LCs. 



than those of the TPC model for the same a and £• (Note, we assume that F^((^, ^>)/i ? 7j tot is 
approximately equal to the ratio of observed energy flux vs. total energy flux.) 

" to 



We can now use the (a,£)-contours obtained by the method discussed in Section 3.1 



constrain the value of fn by computing fu(a, () and overplotting those contours (see Figure 



4. Results 

We infer values for a and £ for each pulsar with typical errors of ~ 5° (see Table [l). The 
tabulated values are the average values of a, £, and (3 implied by the solution contours, while 
the errors are chosen conservatively so as to include the full (non-rectangular) contour. We note 
that once a solution is found at a particular a and £, it is worthwhile to study the LC at the 
position q' = C and (J = a. (i.e., when the angles are interchanged). Such a complementary 
solution may provide a good fit in some cases (see the alternative solutions listed in Table [T]) , 
due to the symmetry of the model as well the symmetry in pulsar geometry (i.e., \(5\ remaining 
constant under this transformation). 

As mentioned in Section[2j we see a shift in best-fit contours, with the TPC model consistently 
yielding smaller a and Q values compared to those obtained using the OG model due to a different 
gap geometry. This is because similar phaseplots are obtained when a larger a is chosen for the 
OG model than for the TPC model. Reproduction of the radio pulse shape then necessitates a 
larger value of £ for the OG fit, since the contour is constrained to appear in one of the grey 
regions (as determined by peak multiplicity; see Figure [I]). The net effect of this is that the 
OG contours are found farther from the origin at a = Q = along the grey bands than the 
corresponding TPC contours. 

Our results for for each of the pulsars (and models) are given in Table [l] 



5. Discussion and Conclusions 

The good constraints on a and £, obtained using the method described above, emphasizes the 
merit of a multiwavelength approach. For example, the light grey regions in Figure [T] indicate 
single-peaked radio profile LC solutions. These would represent typical constraints one would be 
able to derive for a and £ when considering only the radio profile shapes. Similarly, when only 
considering the (single-peaked) 7-ray LCs, one would find relatively large (a,C)-contours [TT] . 
The requirement of fitting both the radio and 7-ray profile shape, as well as their relative phase 
lag, results in much smaller contours, as shown in Figure [T] 

Generally, our best-fit (a,C) compare favourably with those inferred by pQ for the first three 
pulsars (see Table [TJ. However, comparison is hampered by uncertainties in estimating the half 
opening angle, p, which sensitively influences the optimal solutions obtained by [lj. Even a small 
error of 5° on p leads to relatively large errors on the allowed a, so that our best fits would then 
be included in their inferred parameter ranges. Direct comparison of our results with those of 
[I] for PSR J1420-6048 and PSR J1718-3825 is not feasible due to different approaches for 
estimating the radio pulse width, W. We only model the most prominent radio peak, while [1] 
measures W as the width of the total two-peaked profile. Implementing a factor 2 difference in W 
will change the p-contours and thus their best-fit solution considerably, improving the agreement 
between the two approaches. Comparison with [1] is not possible for PSR J1509-5850 as the 
lack of polarization data inhibited inference of a best-fit (a,/3) by [1] in this case. 

The next step for this kind of study is to apply a mathematically rigorous method of 
determining the best-fit LC solutions as a function of many free model parameters to the six 
pulsars modelled in this paper. An example is the Markov chain Monte Carlo method which 
has been successfully applied to millisecond pulsar LCs p2]. 



Table 1. Values for a, £, (3, and fa derived for 6 pulsars for both the OG and TPC models. 
The *'s indicate alternate best-fit solutions upon interchange of a and £. 



Pulsar 


Model 




C(°) 




fa 


J0631+1036 


OG 
TPC 


74 ±5 
71 ±6 


67 ±4 
66 ± 7 


- 6±2 

- 5±3 


0.93 ±0.06 
1.04 ±0.04 


J0659+1414 


OG 

TPC 

TPC* 


59 ±3 
50 ±4 
38 ±1 


48 ±3 
39 ±4 
50 ±4 


-12 ±5 
-13 ±6 
11 ±4 


1.16 ±0.53 
1.64 ±0.04 
1.63 ±0.05 


J0742-2822 


OG 

OG* 

TPC 


86 ±3 
71 ±6 
64 ±8 


71 ±5 
86 ±4 
80 ±4 


-16 ±6 
16 ±6 
15 ±6 


0.99 ±0.10 
0.81 ±0.09 
0.88 ±0.41 


J1420-6048 


OG 

TPC 

TPC* 


67 ±5 
64 ±6 
42 ±5 


45 ±7 
43 ±8 
63 ±5 


-22 ±9 
-21 ±9 
21 ±9 


0.77 ±0.13 
0.90 ±0.10 
0.77 ±0.06 


J1509-5850 


OG 
TPC 


66 ±4 
61 ±5 


50 ± 7 
44 ± 7 


-18 ±8 
-18 ±8 


0.77 ±0.11 
0.89 ±0.10 


J1718-3825 


OG 

TPC 

TPC* 


67 ±6 
61 ±5 
42 ±6 


48 ±6 
43 ±6 
62 ±5 


-19 ±8 
-19 ±8 
19 ±7 


0.76 ±0.12 
0.86 ±0.07 
0.83 ±0.10 
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